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Abstract: The non-perturbative autonomous renormalization of the scalar <E> 4 -model is applied 
in the framework of stochastic quantization. I show that this requires a selective, momentum- 
dependent renormalization of the Onsager coefficient A, a direct consequence of the charac- 
teristic wavefunction renormalization applied. As a result, I obtain a Langevin equation for 
the renormalized constant mode of the field, which is solved numerically. It is demonstrated 
for temperature zero that, starting from specified initial conditions, the system relaxes to its 
equilibrium state, the symmetry-breaking vacuum of the "static" <3? 4 -theory. 
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Recently there has been much interest in the autonomous renormalization (AR) of the Higgs 
field ||] . The AR emerges very naturally as an alternative solution to the renormalization group 
equation [||, 0], where the major difference to the conventional perturbative procedure is an 
infinite wavefunction renormalization |||, [|]. Concerning the number of free parameters the 
AR turns out to be very economical; there is only one for the scale invariant theory . As a 
consequence, it is possible to make a concrete prediction for the mass of the Higgs boson, which 
comes out surprising ly heavy, at about 2 TeV @, |, For more technical and phenomenological 
details concerning the AR, I refer to the literature cited above. 

The question posed in this Letter is whether the AR can also be applied in the context of 
stochastic quantization. The basic idea of this method is to assign a purely relaxational dynamics 
in a fictitious time coordinate, which I will simply call time in the following, to the system [B]. 
Starting from some non-equilibrium initial state, on account of ergodicity and detailed balance, 
the system relaxes to its equilibrium state, the ground state of the euclidean quantum field 
theory. The latter will be called the static theory below. The dynamics is described by a 
Langevin equation and has many formal analogues in statistical physics. There, however, a real 
time-dependent process, like the dynamics of a ferromagnet close to its Curie point, is described 
by the equations. Reasons for introducing stochastic quantization have been advantages in gauge 
theories [|| and numerical simulations Q. 

Consider now the scale invariant euclidean $ 4 theory in a finite periodic hypercube with 
volume Q = L d , described by the action 



(1) 



From the static calculation one knows that the AR requires a selective wavefunction renormaliza- 
tion, i.e., one has to distinguish between constant and finite-momentum (FM) modes Hj. While 
the FM modes are essentially not renormalized, the (bare) constant mode has to be rescaled 
by a logarithmically divergent Z- factor. In a previous work on the static theory [||], it was 
demonstrated that the most natural starting point for the AR is a finite system. Although one 
is not really interested in finite-size effects, many features of the model, like the convexity of the 
effective potential, the emergence of spontaneous symmetry breaking, and finite-temperature 
effects, can best be understood by starting from a finite system and then taking carefully the 
infinite-volume limit. Much the same strategy will be pursued in this Letter for the stochas- 
tically quantized model. For the sake of simplicity, the calculations will be restricted to the 
zero-temperature case. To treat the more general situation, one dimension of the volume f2 has 
to be kept finite ||. 

Following Consoli and Stevenson ||], I denote (in slightly adapted form) the relation between 
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bare and renormalized Fourier amplitudes of the field by 



B (p,t) = Z(p)3 4>{p,t) = [V3e~a <5 Pj0 + (1 - <J p ,o)l 0(p,t) • 



(2) 



Since I will use dimensional regularization in the following, logarithmic divergences are expressed 
by powers of 1/e (with e = 4 — d) from the start [pTCj ] . The vector p in (|2|) has c? components 
with discrete values = 2irki/L, where ki G X. The factor v3 in (Q) is a result of the static 
calculation O. It comes from a finite rescaling of the constant mode (called zq in [0]), which 
guarantees the proper normalization of the (static) propagator. 

The dynamics assigned to the system is described by the Langevin equation 



d t 4>B(p,t) + X B (p)-~ 



5S[(f) 



B\ 



£b(p,<) , 



i>B[-p,t) 

where £b(p, t) is a Gaussian random "force" with mean zero and variance 

(£(p, t) i(Ba, t')) = 2 A fl (p) 5 P ,- q 6(t - t') . 
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From work on classical statistical dynamics [11], it is well known that the renormalization of the 
Onsager coefficient A is closely related to the wavefunction renormalization. Hence, from this 
point of view, it is natural to infer a relation of the form 



A B (p) = [zxf?5 Pt o + (l-6 p<0 )]\ 



(5) 



where the finite normalization constant z\ and the exponent a are to be determined below. 

Like in the static calculation ||, the first step towards a solution of (jH) is the approximative 
integration of the FM-modes in a one-loop type procedure. Similar methods have been applied 
in the context of critical dynamics in finite-size systems 13 1. Here I am following closely the 



method suggested by Goldschmidt p3| . After integrating out the FM modes, one obtains the 
Langevin equation for the constant mode: 



d t Mt) + f (3 c(t) + Mt) 2 ) Mt) = £o(t) 



(6) 



where I have defined 



such that 



bo(t) = - I d d xM^t) 



and £o(t) is related to ^(0,t) analogously. 

In @, the term C(t) represents the effect of the FM modes. It is the time-dependent 
analogue of the tad-pole graph of the static one-loop calculation j|] and may be represented in 
the form 



C(t) = 2XL- 



[ ds^ exp 



-2A 



da' 



t-s 



! i 9 lit l\ 



(7) 
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Suppose that for large time <fio(t) approaches a constant equilibrium value. Then (0) reduces to 
the well-known static result 

C(t = 00) = L- d Y -= -4 ztti — • (8) 

1 ' pToP + 9 Mt = oo)2/2 W 

For the further evaluation, ([?]) may be rewritten as 

C(t) = 2\L~ d J* dsB^^J exp(-\g J° ds'Mt-s') 2 ^) (9) 

with 

00 

£(x) = „4(x) d - 1 and .A(x) = £ e" fc2:r . 

fc=— 00 

The fact that the sum in (0) does not extend over zero momentum has no influence on the UV- 
behavior. Consequently, C(t) contains the typical 1/e pole. In (0) it is caused by the divergence 
of the integrand for small argument, where B(x) ~ (ir/x) d / 2 . 

At this point a brief comment on the initial conditions is appropriate : For the evolution of 
the field expectation, one initial condition will be to assign some specified value to the constant 
mode at t = 0. Moreover, from (0) and (p|), it becomes obvious that another "initial condition" 
is C(t = 0) = 0, i.e., fluctuations of the FM modes are suppressed at t = 0, which, in a sense, 
emerges here as the natural initial state of the relaxational process. 

For the further evaluation of (^) it is convenient to isolate the pole by rewriting the whole 
expression in the form 

C{t) = Cdiv(t) + AC(t) + O (L- d ) . (10) 

The contributions O (-^ do not influence the dynamics in the infinite-volume limit. 
The UV-divergent part Cdi v is given by 

C div(t) = { *>/2 I dss ^ /2ex P (-A5^o(i) 2 ) • (11) 

The integral can be calculated exactly, and, after analytic continuation to 4 — e dimensions, the 
result may be expanded in a power series in e : 

(4vr) 2 e 

where the finite contribution is given by 



C d Ut) = - y -^ + C2(t) + 0(e), (12) 



r (°) - e e ■ 9 4>l 



(4vr) 2 2At 2(4vr) 2 



log Vlo 8 (6)e-e d ^' 8 > 



16vr 2 $ 2 / fa v ; dx 



x=l 



(13) 



with the abbreviation = @(t) = Xtg4>o(t). In (|l3|), j(x,y) denotes the incomplete Gamma 
function and Ce Euler's constant. Q v emerges in this stage of the calculation as an arbitrary 
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momentum scale generated by dimensional transmutation, but its normalization is chosen such 
that it later on also turns out as the vacuum expectation value of the field. 
The second term on the rhs of ( |i~0| ) is a memory term given by 

a-d/2 f t 



(14) 



AC(t) = M_ jT* dss- d / 2 [exp (-XgJ^ds'Ms'f) -exp (-XsgMtf 

Like in the static calculation ||, the pole in (|ll]) is cancelled by the nonlinear term in (^) 
by assigning the UV-flows 

Mt) = VSe- 1 / 2 ^) and g = (15) 

to the constant mode (cf. Eqn. (0)) and the bare coupling, respectively. After this cancellation, 
the remainder cj^(t) + AC(t) on the lhs of @ is of order e°. Thus, together with the factor 
g, the deterministic "force" would be of order e. A finite contribution may be obtained by 
choosing a = —1 in (||). Moreover, the normalization z\ can be fixed by demanding that the 
deterministic part is equal to the derivative of the static effective potential for t — > oo. In this 
limit, the memory term AC(t) tends to zero, and only the first term in square brackets in (|l3| ) 
survives. As a result, one derives the relation 

from which it is also obvious that = $> v (where the logarithm has its zero) corresponds to the 
asymmetric ground state. Now the second derivative of V with respect to the renormalized field 
must be equal to the squared propagator mass, called Mh in the following. The latter can be 



read from (g), and inserting the UV-flows (15) one finds 



M 2 U = 8vr 2 §1 . 

Eventually, the second derivative of V, calculated from (|i~6|), is equal to Mfj if z\ = 3. 

What, then, happens to the noise in view of finite and infinite renormalizations of the field? 
Carrying out the rescalings in (||) gives rise to a renormalized Langevin equation of the form 
dt&(t) + finite = 3 _1 / 2 e 1 / 2 ^(t). Thus, the most natural UV-flow for the constant mode of the 
noise is 

The renormalized noise S(i) has mean zero, and with (0) and (||) one obtains for the variance : 

(E(t)E(t')} = 2\L- d 5{t-t'). (18) 
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The average magnitude of the noise is 0(L~ d l 2 ), and, therefore, in the infinite- volume limit the 
dynamics becomes essentially deterministic. This is completely analogous to the static case, 
where a saddle-point approximation for the path integral over the constant mode becomes exact 
in the infinite- volume limit |J. 

With (|l3|), djjj), and ( |l5| ) in (|6|) and after introducing dimensionless variables i 
and = &/& v , the equation of motion for the constant mode reads 

d7 (x,i$ 2 



2XtM 2 H 



log(t^)exp(-t<^) + 



dx 



+ C E - log($' 



x=l 



+ 



exp(— t $ 2 ) 



ds s 



exp 



ds'$(s') 2 -exp 



(19) 



where the argument of $ has been suppressed when it is t. 

There is one special feature, which one encounters on the way to a numerical solution of 



(19). The first term on the rhs of (|i3| ) and, in turn, the first term in the second line of ( |l9f) are 
singular for i — > 0. This is due to the initial state discussed above, and the singularity could 
be avoided by choosing different, though more complicated, initial conditions. Similar so-called 
short-time singularities are well known from statistical physics, where they give rise to anomalous 
behavior of relaxational processes in their early stages [jl4| . In stochastic quantization, where 
one is mainly interested in the equilibrium state, the short-time behavior seems to be merely a 
technical problem. A non-singular equation may be obtained with the ansatz = i 1//4 ^(t ) 
with finite initial value ^(O). 

A number of numerical solutions is depicted in the Figure. ^(0) serves as a parameter. As a 
consequence of the short-time singularity, all trajectories of <I> have to start from zero. Then the 
expectation value increases and overshoots the equilibrium value, but asymptotically, for large 
time, all solutions approach = 1 from above. There are no solutions that cross the origin. 

In summary, I have applied the autonomous renormalization to the stochastically quantized 
<& 4 -theory. The most important difference compared with conventional perturbation theory is a 
selective renormalization of the Onsager coefficient A as formulated in Eqn. (||). A finite equation 
of motion, which becomes purely deterministic in the infinite- volume limit, has been derived and 
solved numerically. Asymptotically, all solutions approach the equilibrium state, the asymmetric 
ground state of the static <I> 4 -model. 
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Figure captions 

Fig. Numerical solutions of Eqn. ( |l9|) with the ansatz = f 1//4 ^ / (t), where ^(0) serves as 
the parameter. 
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This figure "figl-l.png" is available in "png" format from: 
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